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A METHOD FOR THE CALCULATION OF EXTERNAL LIFT, MOMENT. 
AND PRESSURE DRAG- OF SLENDER OPEN-NOSE BODIES OF 
REVOLUTION AT SUPERSONIC SPEEDS 
By Clinton E. Brown and. Hermon M . Parker 

SUMMARY 


An approximate method is presented for the calcula- 
tion of the external lift’, moment,, and pressure drag of 
slender open-nose bodies of revolution at supersonic 
speeds, The lift, moment, and pressure drag of a typical 
ran- jet body shape are calculated at Mach numbers of 
1,1x5, l.oO, I. 75 , and 5,0 and the. lift and moment results 
are compared with available experimental data. The. 
agreement of the calculated lift and moment data with 
the experimental data is excellent. The pressure -drag 
comparison was not presented because of the uncertainty 
of the amount of skin-friction drag present in the 
experimental results. It was found that the- lift coeffi- 
cient definitely increased with increasing Mach number,.- 
whereas the moment coefficient taken about the midpoint 
of the body and the drag coefficient decreased with 
increasing Mach number. The manner in which : the method m 


annular air inlets is Shown, 


be reliably used for obtaining the aerodynamic charac- 
teristics of slender bodies that are required for effi- 
cient supersonic flight. 


es 

of 

revolut j 

Lon 

v/i 

th 

The 

e' 

xcellent 

agree: 

merit 

eht 

r 

e suits w:i 

Lth 

ex 

peri- 

ap 

>pr 

oxirnate v. 

letr 

lod 

may 

, i-'-. 

; 0 -: . 

e 

aerodynai-i 

lie 

eh 

arac- 


INTRODUCTION 


Current proposals for the design of aircraft capable 
of sustained flight at supersonic speeds and utilising 
the ram jet as a method of propulsion nave established 
the importance of 'knowing the aerodynamic characteristics 
of slender one-i-nose bodies of revolution at 


speeds 
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greater than that of sound. The lack of theoretical 
treatments and experimental data, emphasizes the need for 
theoretical investigation of this problem to serve as a 
guide for future work and as a check oh the reasonable- 
ness of current and future experimental results. 

The small "perturbation approximation was used in 
reference 1 to deduce the wave drag and in reference 2 
to obtain the lift and moment of slender pointed-nose 
bodies of revolution. No fundamental analysis is known 
to have been made, however, of the characteristics of a 
slender open-nose body shape, such as that required by 
ram-jet propelled craft. The peculiarity of the problem, 
from general considerations of similarity, is that the 
flow pattern is two-dimensional at the lip of the nose 
and approaches the three-dimensional pattern farther 
along the body. The present work extends the method of 
references 1 and 2 to apply to these slender open-nose 
bodies of revolution with supersonic flow into the nose. 

'The result is a fairly simple method of numerical inte- 
gration of the differential equation of the flow. As an 
illustration, the pressure distribution, wave drag, lift, 
and. moment are calculated at Mach numbers of l,i[„ 5 > I.60, 

I.75, and 5.0 for a typical ram- jet airplane body shape 
and the lift and moment results are compared with the 
experimental data. It should be pointed out that the 
accuracy of the method, which assumes potential supersonic 
flow throughout the field and also assumes small disturbances, 
depends on the surface angles of the body and the Mach 
numbers The error increases with either increasing Mach 
number or increasing surface angles. 


SYMBOLS 


x, r, 0 
X 

l 

R 


P 


cylindri cal coordinates 

distance along X-axis measured from nose of 
body 

length of body 
radius .of body 
Mach angle 



b = / a 


1 
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j/(x, v f 0) 

perturbation potential 

(x, r) 

perturbation potential for axial flow 

^ 2 (x, r, 0 ) 

perturbation potential for cross flow 

v x 

/6cf\ 

axial velocity increment \&J 

v r 

radial velocity increment 

V 

velocity in undisturbed stream 

a 

velocity of sound in undisturbed stream. 

M- 

Mach number in undisturbed stream (v/a) 

P 

density in undisturbed stream 

Ap 

incremental surface pressure due to angl 
of attack 

Pi 

local pressure 

;o 

pressure in undisturbed stream 

Y 

ratio of specific heats of air (ijp) 

a 

angle of attack, radians (except where 
o the rw i s e no t e d ) 

6 

angle between surface of body and X-axis 

°L 

lift coefficient ^Lif t TTR-^^ . 

Cp 

drag coefficient (Drag^^'trR^J 

C m 

moment coefficient ^Moment 

u 

variable of integration 


X 


Br cosh u 


J'i - x j _ - BHj 
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n 


n 


BR. 


n 


f '<£). 


Subscr>ipt3 and superscripts: 

H refers to nose 

n refers to n^ n integration station, summation 

variable 

i refers to i^ n integration station, summation 

variable 

deg in degrees 


MATHEMATICAL .ANALYSIS 
Pointed Bodies of Revolution 


The line ariaed equation of motion of a nonviscous 
compressible fluid may be written for a cylindrical 
coordinate system: 


dr 2 r Or r 2 dQ* N ' bxr 


( 1 ) 


■where $ is the potential function assumed to represent 


the effect of a small disturbance set up by the slender 
bodies being considered. The problem is to find a solu- 
tion of equation (1) that will satisfy the known boundary 
conditions at the surface of the body. A general solu- 
tion of the differential equation (1), when M > 1, for 
diverging waves has been found by Lamb (reference 3) to 
be, with a slight change in notation,. 


0 - X - Q s r s cos 


P s r s sin s0 


where 
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Of? — 


ana 


% 


*o 


f\CO 

0 

f \ oo 


m 




( 2 ) 




0 


y 


f (x - Br cosh u) du 


0 


g(x - Br cosh u) du 


where 


B — Vk 2 “ 1 


The part of Lamb * s general solution corresponding to con- 
verging waves does not apply to the present problem 
because all disturbances originate on the body and diverge 
into the flow field investigated, yon Karmah and Moore 
have investigated the problem of the resistance of pro- 
jectiles and cones (reference 1} and have found a solu- 
tion for the case of axial symmetry 


n» 


0 


f(x - Br cosh u) du 


( 3 ) 


o 


a 

In 


case of the general 
analysis the body, in 
vojectilo , was represented by 

starting at 


sped a] 
their 


a 


that can be seen to be 
solution with s = 0. 
this case a sharp-nose 

distribution of sources along the X-axis 
x ss 0, the nose of the body. By a numerical method of 
integration, it became possible to write the equations 


for the 


velocity 


increments v-p 


and 

& 


V 




V, 


’ll 


n 

-*Y 




L=1 


T - 1 \/( T i-i ri ) • ~ 1 


( 4 ) 
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n 


v 


x n 


1=1 


= X cosh" 1 (T i n ) - cosl-r 1 (r 1 _ 1 n ') 


( 5 ) 


where 


x„ - J. 


and 


m n n J r 
A i. “ ~ C o 
pil n 


3i X 1 ~ 3^1 


A ± = f » (&) i 


with the boundary conditions 


v. 


V + v x 


dr 

dx 


( 6 ). 


These three equations , in three unknowns v r n ? v Zn) 

were solved at each station on the body for v r and 'v 
The pressures were then found from the Bernoulli equa- 
tion in the ' fo.ru: 


p 




-X . 

*1 

Y -x 


( 7 ) 


Ferrari (reference Lj.) and Tsien (reference 2) have 
independently found solutions for the case of pointed 
bodies of revolution at small angles of attack. Their 
solutions showed that the potential could be expressed 
in two terns : the firs^, from equation (5 ) , 

pco 

~ - fq(x - Br cosh u) ' du 

p 0 
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is the solution for the pure axial flow already described, 
and the second 


- -3 cos 8 f^(x - Br cosh u) cosh u du (8) 

Jo 

represents the cross-flow potential of an arbitrary dis- 
tribution of doublets along the axis of the body starting 
at the nose of the cone- or projectile. The form of 
equation (3) is such that the cross flow is from the 
direction 0 5 0, as shown in figure 1. 

By neglecting the small effect ox* the axial flow on 
the lifting pressures, Tsien obtained for the pointed 
projectile of arbitrary shape the equations: 



•^a+I ~ x n 
R base 


n 


R n+X +R n\ 


R 


base 


~L 


. i«0 


2 

III 

2 clY 



( 9 ) 


c r , — — ' 


2a\ 

x n-t~l + 2x n 

m / 

3 

n 


x 1 1 - 
n+1 

x r„ ■+ : 

n u+l 


x. 


+ 


n+I ~ vV n 


R n+1 


R + R 
n+1 h , 


n 


R 


base 


R 


base 




2 ay 
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r> 

C 


n 


1 


i-l 


■p2xr « 
S3 j. >. “j 


2 a? 



- cosh 



J C 




( 11 ) 


The values Kj_ In these equations are assumed to be 
constants for each' .interval of the. step-by-step process . 
The moment coefficient of equation (10) is assumed posi- 
tive for nosing-up moments, these moments being taken 
about the nose. 


Open-Nose Bodies 

The flow conditions over an open-nose body, differ 
from those of pointed bodies in that, for finite angles 
of the nose lip, the flow is two -dimensional at the' lip. 
This problem was not considered in references 1, 2, and 1}. 
and the general solution should therefore be examined to 
determine its applicability to this special case. Lamb 
has shown (reference 3) that a sufficient requirement 
for the existence of the general solution to the dif- 
ferential equation of motion .{reference 1) is that 
f(x - Br cosh u) be zero "for all values of the argu- 
ment less than some arbitrary limiting value. Hie deter- 
mination of f(x - Bn cosh u) such that the boundary 
conditions at the open-nose body are satisfied assures 
the fulfillment of this general requirement* .For the 
usual case of supersonic flow into the nose, the boundary 
condition requires the surface of the body to be a con- 
tinuation of a cylindrical stream surface of radius R- k j 
in the undisturbed flow ahead of the body as shown in 
figure 2. The perturbation potentials, equations 0) 
and (8), therefore must be zero at the cylindrical 
stream surface ahead of the body. Substituting 
£ - x - Br cosh u in equations' (3) and (8) gives 


-Br 




X q 


0 


zr 


B?r 2 


— 00 


(12) 
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and 


p x-Br 


$ 2 . 


cos 0 


W — oo 


f 2 (§)(x - £)d£ 
\/(x - £} 2 ~ 3 2 r 2 


U3) 


It Is obvious 
by letting 


f 1 ^ ^ 


the boundary conditions are satisfied 
: f P (c, ) = 0. for all values of 


g < Xq - BRq, where the point 


the open-nose 


y< 


( x o> R o) 


at the lip oi 


It then remains to determine f-, (§) 


and f 2 (& ) for S, > Xq - BRq so that the body surface 

Is a continuation of this stream surface. From physical 
considerations, f^ (.£,) and fp(£.) may bo regarded as 


an axial distribution of sources and doublets, respec- 
tively, where s is measured along the X-axis. Because 
the effect of a source or doublet can be felt' only along 
or behind its Mach cone, the., source distribution must 
begin a distance BRq ahead of the nose. This point is 

chosen for the origin of the coordinate system. (See 
fig. 1. ) It must be emphasized that the source and 
doublet distribution determined by satisfying the 
boundary conditions at the stream and body surfaces shown 
in figure 2 does not represent correctly the flow inside 
that stream surface. This result corresponds to, the 
physical fact that the actual supersonic flow into the 
nose does not affect the flow external to the body. The 
basic assumptions of potential flow and small disturb- 
ances are valid provided the slope of the body surface 
is small. Actually., for finite angles of the nose lip, 
a nonconical shock wave is formed that causes a loss in 
total head and produces rotation in the field. 


Numerical integration of equation (12), with con- 
stant values of f-, 1 (£,) assumed over the Integration 

intervals, results in the same expressions for Vv, 

• r n 

and v-,, as those obtained by von K arm ah and Moore 
x n • 

(equations ( Ip ) and (5))« -These constant values of 
f •(£) are determined by satisfying the boundary 

condition; 
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v r 

= tan 5 (lit.) 

V + v 11 

X 

where tan 5,_ is the slope of surface of the body at 
the n u ‘- interval of integration. 

By following the method of reference 2, the lift and 
moment coefficients for* small angles of attack based on 
the area of the nose may be written 


PIT P 


C L 


1%V 


0 o 


3P, 


N '6^2 « ,a , 

'T» fca - COS: 0 do R ax 
ox 


H 


nTT 


Cm" 




TTRij^V „ 


p l+BR 


0 


iM 




{1 \ U F; 

(x - g* - BR^y cos 9 d0 R dx (16) 


ER,r 




where l is the length of the body, R^ is the nose 
radius, and' the moments are taken about the midpoint of 
the body. By substituting the expression for dj^^/dx 

in -equations. ( 1R ) and (1.6) , CL and b-ec-ome 


r\ TT 


C E 


•m-B 


•HR 


n 2 V 


2 q 
cos 0 


PZ.+3R 


n 


de 


R dx 


u 


0 


J BR 


N 


oX^BR 

f 2 '(§)(*■- g)dg 
0 W(X £) 2 ~ 3 2 R 2 



ai 
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Cra = 


4 


VJ J 


r>TT pl+BR- 

. cos 2 0 . d9 

0 ' °B% 


H 


(x - | - Br)h dx 




•BH 


fV(s)U - £)d£ 


0 


BP 


l-ycT-.S) 2 - B 2 R 2 


(i8) 


The distribution function f 2 (S ) must be determined by 
the boundary condition; thus, 

pa cos 9 = ( ~^) (19) 

\ 6r / p—p 


for which the radial velocity is assumed to be normal' to 
the surface', A more rigorous boundary condition taking 
into account the slope of the body was giver, by Ferrari 
(reference 1|). For small surface angles, however, 
equation ( 19 ) is within the accuracy of the small- 
perturbation assumptions. The expression 


( 



cos c 


R 


:-BR 


^0 


f 2 »(g,)(x - £) 2 dg 
'JiX - l) 2 - 3 2 R 2 , 


( 20 ) 


is integrated numerically for constant values of 
f 2 1 ( § ) = over the i tn interval of integration to 

obtain the sum 
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Substituting this equation in equation (19) gives 


n 



With the values of - — determined, equations (17) 

2Va 

and (l8) become 
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In equations (23) and (2Lj.) the pressure used for a 
given integration interval is the average of the pres- 
sure 3 at the beginning and at the end of the interval'. 
This scheme of using average lifting pressure is par- 
ticularly necessary in regions where the pressure is 
rapidly changing. The method does not give the pressure 
at the beginning of the first integration interval, that 
is, at the point h - 0, It can be shown that, as the 
first interval approaches zero, the pressure at the 
lip (n s= 0) is obtained by letting the expression in 
equations ( 23 ) and (2l>.) 



i=l 

have the value 0,5 when n ~ 1. 


Method of Calculation 

If calculations are to be made for an open-nose, 
body, the total number of integration stations chosen 
is a compromise between the amount of labor involved and 
the accuracy desired, due consideration being given to 
the limitation on the accuracy imposed by the basic 
assumptions of small disturbances and potential flow. 

In general, where the pressures are changing most rapidly 
the integration stations should be the most dense. 

The integration stations must first be chosen. (See 
fig. 3.) The calculation for .the lift and moment then 
proceeds as follows : 

The boundary condition, equation (22), is applied 
to the point n = 1, the summation reducing to a single 
term, where 




j 


1 


BR 


1 
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and 


T 


n 

i-I 


H - 0 
BR 1 


Equation (22) 

equation (22) 
now two terms 


then gives a value of B^Kq/^Vu. Next, 

is applied to point '2, n = 2, . containing 
Substitution of the expressions 


T?. 


A 2 


- J 2 


BR: 



= x 2 " J'l 
br 2 


B 2 K 1 /2Va permits 
the calculation of the value of B^K 2 /2V a. 

Bp continuing the process, at each successive 
station one more term occurs in equation (22) involving 
one new K-?., which is then determined. With the values 

of 3 K^/2va determined, . equations (25) and (2 Ip) are 

used tc evaluate the lift and moment coefficients. It 
will be noticed that, when equations (25) and (21*.) ‘are 
evaluated for n = 1, the expression occurring after the 
last summation symbol is unobtainable. For reasons 
previously stated, the expression must be given the 
value O.'y. 

The procedure f or •calculating . the drag pressures' is 
similar to but somewhat simpler than that for -the lift 
and moment. Equations (l;.) and (5 ) are evaluated at the 
point n - 1, the sums reducing to one terra. By applying 
the boundary conditions (equation (6)) and using the 


m 2 - 12 
"0 


- 0 


3R 2 


and use of the value obtained for 
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known slope of the body dr/dx, the constant A-^ occur- 
ring in equations (I4.) and ( 5 ) is determined. Substi- 
tution of A] back into equations (I 4 .) and ( 5 ) gives the 


increment velocities at n = 1. Bernoulli’s expression 
(equation (?)) then gives the pressure ratio at n - 1> 

It is to be noted, since V x Ma, that, actually, values 


of Aq/a , 


v. 


r l 


j a. 


and /'a are determined. Pro- 


ceeding to the point n ~ 2, one new A is involved in 
equations (It.) and ( 5 ) that is determined by the boundary 
condition in equation (6). In the same way, the 
velocities and the pressure ratio are calculated at 
n ~ 2. The process is continued at successive integra- 
tion stations over the whole body, The calculated values 
of the pressure ratio allow a pres sure -distribution curve 
to be drawn and the drag coefficient to be evaluated in 
the usual manner, 


DISCUSSION OP RESULTS 
Results of Calculations 


Calculations were made in order to obtain the pres- 
sure drag, lift, and moment of a typical open-nose body. 

A sketch showing the dimensions, the integration stations, 
and intervals is given in figure 2 . Calculated pressure 
distributions at zero angle of attack for the Mach 
numbers l,k% I. 60 , 1 . 75 * an & 3 #00 are presented in 
figure k. The pressure rise at the nose lip is approximately 
that which would be obtained over a two-dimensional wedge 
of the same angle. Ackeret 1 s theory for small, disturbances 
(reference 5 ) gives for the pressure rise of a 5 0 wedge 

p 2 ” 

— » l.lkT at M ~ l.k5, whereas the pressure ratio on 

P 1 

the lip of the open-nose body is about 1 . I 4 O, This 
agreement is considered a reasonable check inasmuch as 
the pressure on the lip must be an extrapolation of the 
pressure distribution from the first point of the integration 
process, which must _be a finite distance back of the lip 
edge. It can be shown that, as the size of the first 
interval approaches zero, the pressure at the lip becomes 
that given by the two-dimensional theory. The effect 
over the nose section of the size of integration inter- 
vals is illustrated by figure 5 * The lYrpoint method 
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was chosen because a greater number of points resulted 
in only a small increase in accuracy and a large increase 
in the labor involved. As would be expected, for the 

case of the. straight conical nose, the pressure falls off 
from the leading edge and approaches the pressure of a 
cone of the same surface angle (reference 1). At. the ■ 
corners of the body, the pressure falls approximately in 
accordance with; the Prandti- Meyer relation (reference 6) 
for supersonic flow around a two-dimensional corner. On 
the center and tail sections there are positive pressure 
gradients that for actual flight conditions would tend to 
cause separation. The- source distribution functions 
corresponding to the ores sure distributions of figure ij. 
have been plotted in figure 6. 


The incremental surface pressures giving rise to 
the lift and moment- are shown in figure 7 lor the Mach 
numbers 1.6 and 3 . 0 . At the higher Mach numbers the 
pressures decrease less rapidly over the nose section. • 

The doublet distribution functions at the four Mach num- 
bers are presented in figure 8,’ The curves of figure 9 
show the theoretical variation of lift-curve slope, 
moment-curve slope, and drag coefficient with Mach number. ’ 
The rather' interesting result obtained is that the lift, 
coefficient increases with Mach number * Experimental 
results presented, in reference 7 show an increasing lift- 
curve slope with Mach number for a pointed projectile. 

The fact that the .doublet distribution for the open-nose 
body is similar to that of a typical projectile (refer- 
ence 2) leads to the expectation that the Mach number 
characteristics of the two shapes would be similar. The. 
center' of pressure at the lower Mach numbers is ahead of - 
the nose and moves back with increasing Mach number. The 
drag coefficient can be seen to decrease with Mach number 

1 

but to a lesser degree than the ~ law; of twro- 


dimensional wing profiles. 


Comparison with Experiment 

A comparison of. the calculated Lift and moment coef- 
ficients with some experimental data obtained in the 
Langley 9~i n ch supersonic wind tunnel is presented in 
figures 10(a), 10(b), and 10(c). The contribution to the 
lift of the internal air can be showm to be; 

ACl = 2a 
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This increment will appear at the nose and will there- 
fore give rise to a 'moments 

AC m = a 

Hies e increments have been added to the calculated 
results for comparison with the experimental data that 
include the effect of the internal air. The experi- 
mental lift and moment results did not go through the 
origin because of deviations in the stream flow and 
perhaps small, unaccounted for, tare forces. It is 
obvious that for bodies of revolution the data should 
give symmetrical curves about the origin. The ’ experimental 
results were therefore made to go through the origin by 
subtracting from each point the value of the lift coeffi- 
cient or moment coefficient at zero angle of attack. 

All coefficients are based on the nose area. 

For it ~ 1 .It 5 two calculated curves corresponding 
to two body shapes tested are shown. The two shapes, 
which are designated for convenience regular tail and 
ferrule tail, differ only in' that a stabilizer-supporting, 
ferrule was placed over part of the tail section. The 
two configurations are shown in figure 10(a). No dif- 
ference can be observed in the experimental results for 
the two shapes, whereas the calculations show a small 
difference. A possible conclusion is that thickening of 
the boundary layer ahead of the shock at the trailing 
edge tends to make the effective shape of the regular 
body more nearly like that of the body with the tail 
ferrule. The better agreement of the ferrule tail cal- 
culation is further evidence that this thickening of the 
boundary layer actually occurs. In any case, the cal- 
culated results are in excellent agreement with the 
experimental data. At M = 1. 60 and M .= 1.75? only 
the regular shape was calculated, and, although the 
agreement with the experimental data is good, ferrule 
tail calculations would probably have given even better 
agreement. The discrepancy at M~ 1.75 between the 
theoretical and experimental points near ±5° ts probably 
due to, separation phenomena. The agreement is good at 
small angles of attack. Because of the uncertain value 
of the skin-friction drag, present in the experimental 
data, a drag comparison is not presented. 
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Another Application 

The success of the method when it is applied to 
open-nose bodies suggests its use in the ' calculation of 
external lift, moment, and pressure drag of bodies of 
revolution having annular air inlets (fig. 11 ) . The 
characteristics of a body of this type would be calcu- 
lated 'as follows: The source and doublet distribution 

from O' to .B would be numerically determined, by the 
methods of references 1 and 2, which would allow.' the 
calculation of the lift, moment, and drag contributions 
of the portion of the body from 0 to C f In order to 
calculate the lift, moment, and drag of tire ..remaining • 
external shell, the source and ^doublet distribution 
from 0 to A would be retained 'and a "new source and 
doublet . distribution .beginning at A and satisfying the 
boundary condition at the outer shell would be ’deter- 
mined. ‘ In this way, the. Influence of the forward portion 
of the body on the flow along the outer shell would be 
fully taken into account. The flow at the .annular Inlet , 
could be fully determined from the original source and 
doublet distribution from 0 to B, which would allow 
the determination of the lift due to the internal flow. 
The moment would be closely approximated by assuming 
this incremental lift to act at the lip of the Inlet. 


CONCLUSIONS 


An approximate method was presented for the calcu- 
lation of the external lift, moment, and pressure drag 
of slender open-nose bodies of revolution at supersonic 
speeds. The pre s sure ^drag comparison was not presented 
because of the uncertainty of the amount ,of skin-friction 
drag present in experimental results. The calculated 
lif t and moment results at Mach numbers of l.ij.5, 1.60, 
1,75, a nd 3,0 showed excellent agreement with the avail- 
able experimental data. The excellent agreement of the 
calculated lift and moment results with experimental 
data indicated that the approximate method majr be reliably 
used for obtaining the aerodynamic characteristics of 
slender bodies that are required for efficient supersonic 
flight. 

Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, ya. 
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Figure 1,- Coordinate system. 
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Figure 4.- Calculated pressure distributions for 
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Fig. 7 



Figure 7.- Distribution of incremental surface pressures 
over the body at two Mach numbers. 
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Figure 9.- Theoretical variation of lift, moment, and drag 

coefficients with Mach number. 
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Figure 10.- Comparison of calculated lift and moment coefficients 

with experimental data. 
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